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ON  THE  SECONDARY  INSTABILITY  OF  THE  MOST 
DANGEROUS  GORTLER  VORTEX. 


S.  R.  Otto 1  and  James  P.  Denier 2 
^CASE,  NASA  Langley  Research  Center,  VA  23681-0001,  USA 
2University  of  New  South  Wales.  Kensington,  NSW  2033,  Australia 

ABSTRACT 

Recent  studies  have  demonstrated  the  most  unstable  Gortler  vortex  mode  is  found 
in  flows,  both  two  and  three-dimensional,  with  regions  of  (moderately)  large  body 
curvature  and  these  modes  reside  within  a  thin  layer  situated  at  the  base  of  the  con¬ 
ventional  boundary  layer.  Further  work  concerning  the  nonlinear  development  of  the 
most  dangerous  mode  demonstrates  that  the  flow  results  in  a  self  induced  flow  re¬ 
versal.  However,  prior  to  the  point  at  which  flow  reversal  is  encountered  the  total 
streamwise  velocity  profile  is  found  to  be  highly  inflectional  in  nature.  Previous  work 
then  suggests  that  the  nonlinear  vortex  state  will  become  unstable  to  secondary,  in- 
viscid,  Rayleigh  wave  instabilities  prior  to  the  point  of  flow  reversal.  Our  concern  is 
with  the  secondary  instability  of  the  nonlinear  vortex  states,  which  result  from  the 
streamwise  evolution  of  the  most  unstable  Gortler  vortex  mode,  with  the  aim  of  deter¬ 
mining  whether  such  modes  can  induce  a  transition  to  a  fully  turbulent  state  before 
separation  is  encountered. 


Research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  contract 

No.  NASI -19480  while  the  author  was  in  residence  at  the  institute  for  Computer  Applications  in  Science 

and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23881 
2 

Research  was  supported  by  the  Australian  Research  Council 


iii 


1.  Introduction 


The  most  dangerous  Gortler  vortex  mode  has  recently  been  identified  by  Denier  et 
al  (1991)  (see  also  Timoshin  (1990)).  This  mode  is  found  to  have  a  streamwise  growth 
rate  of  0(G 3/5)  and  is  localized  in  a  thin,  0(G-1/5),  viscous  layer  located  at  the  solid 
boundary.  The  effect  of  crossflow  on  this  mode  has  been  demonstrated  by  Bassom  Sc 
Hall  (1991)  to  have  a  stabilizing  effect.  The  nonlinear  evolution  of  the  most  dangerous 
mode  has  been  considered  in  a  series  of  articles;  for  two  dimensional  flows  by  Denier  Sc 
Hall  (1993)  and  for  three  dimensional  flows  by  Otto  &  Bassom  (1993).  In  both  cases 
the  vortex  induced  mean  flow  results  in  a  region  of  reversed  flow  at  some  finite  distance 
downstream  of  the  position  at  which  the  perturbation  is  first  introduced.  However, 
prior  to  such  a  flow  reversal  the  total  streamwise  velocity  fields  are  seen  to  be  strongly 
inflectional  in  nature  and  thus  we  anticipate  that  such  flows  will  be  highly  susceptible 
to  secondary  inviscid  instabilities  in  the  form  of  Rayleigh  waves. 

In  the  case  of  order  one  wavenumber /Gortler  number  regime  Hall  Sc  Horseman 
(1991)  found  that  the  highly  inflectional  velocity  profiles  found  in  Hall  (1988)  are 
unstable  to  inviscid  Rayleigh  instabilities  and  demonstrated  good  agreement  with  the 
experimental  results  of  Swearingen  Sc  Blackwelder  (1987). 

In  the  situation  under  consideration  the  equation  governing  the  perturbation 
quantities  is  the  three-dimensional  Rayleigh  pressure  equation  obtained  by  Hall  Sc 
Horseman  (1991).  We  obtain  solutions  of  this  boundary  value  problem  in  order  to 
ascertain  the  temporal  stability  of  a  given  vortex  flow  profile.  Since  the  presence  of 
the  vortex  implies  that  the  underlying  basic  state  has  a  periodic  spanwise  structure, 
we  elect  to  retain  the  spanwise  variation  by  only  assuming  that  the  periodicity  of  the 
secondary  modes  are  the  same  as  the  that  of  the  vortices.  In  fact  we  find  that  the 
most  unstable  modes  are  comprised  of  significant  contributions  from  the  mean,  fun¬ 
damental  and  second  harmonic  components.  Similar  observations  have  been  made  by 
Balachandar,  Streett  Sc  Malik  (1992),  in  their  work  concerning  the  secondary  insta¬ 
bility  of  rotating  disk  flows.  Somewhat  surprisingly  this  procedure  does  not  generate 
any  odd  modes,  those  which  are  t/2  out-of-phase  with  the  underlying  vortex  state, 
of  the  form  found  by  Hall  Sc  Horseman  (1991).  However,  by  explicitly  apply  spanwise 
boundary  conditions  applicable  to  the  odd  modes  additional  unstable  modes  can  read¬ 
ily  be  obtained  but  which  have  smaller  growth  rates  than  the  those  labeled  secondary 
in  Figure  1,  (see  Otto  Sc  Denier  (1993)  for  further  details). 

The  remainder  of  this  article  is  structured  as  follows:  in  section  2  we  shall  sum¬ 
marize  the  derivation  of  the  equations  governing  the  nonlinear  vortex  state,  and  derive 
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the  inviscid  Rayleigh  pressure  equation.  In  section  3  some  brief  comments  will  be  made 
concerning  the  numerical  methods  used  to  solve  the  Rayleigh  equation,  and  in  section 
4  we  shall  comment  on  our  findings  and  discuss  possible  future  topics  of  interest,  and 
finally  in  section  5  we  shall  draw  some  conclusions. 

2.  Governing  equations 

In  this  article  we  consider  a  boundary  layer  flowing  over  a  yawed  cylinder,  this 
flow  has  been  previously  considered  in  Hall  (1985)  and  the  reader  is  referred  to  that 
paper  for  details.  The  Reynolds  number  Re  and  Gbrtler  number  G  are  defined  by 
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where  V  is  a  typical  flow  velocity  in  the  streamwi.se  direction,  L  is  a  characteristic 
streamwise  lengthscale  and  v  is  the  kinematic  viscosity  of  the  fluid.  The  curvature  of 
the  cylinder  is  taken  to  be  (f)  where  the  function  \o  is  supposed  to  be  smooth 
and  positive.  With  these  definitions  S  =  Ljb  where  b  is  a  typical  radius  of  curvature 
of  the  cylinder.  The  Reynolds  number  is  taken  to  be  large  whilst  6  is  sufficiently  small 
so  that  as  S  — *  0  the  parameter  G  is  fixed  and  is  of  order  one  (when  compared  to 
the  Reynolds  number).  We  will  subsequently  consider  the  large  Gbrtler  number  limit 
relevant  to  the  most  dangerous  Gbrtler  vortex  mode. 

The  nonlinear  evolution  of  the  most  dangerous  Gbrtler  mode  occurs  over  an 
((2-3/5)  Spatia,l  lengthscale,  (where  now  we  are  assuming  G  >  1),  and  is  confined  to 
an  0(G-1/5)  viscous  layer  located  at  the  solid  boundary.  With  the  effect  of  crossflow 

—  —  3  2 

of  order  0(Re  2  G * )  the  modes  evolve  over  an  0(G~  * )  temporal  scale.  The  resulting 
equations  governing  the  nonlinear  evolution  of  the  most  dangerous  Gbrtler  mode  are 
given  in  Otto  &  Bassom  (1993)  (also  Denier  &  Hall  (1993),  Timoshin  (1990)).  For  the 
sake  of  brevity  the  reader  is  referred  to  the  aforementioned  papers  for  full  details. 

Here  we  make  a  few  brief  remarks  concerning  the  results  of  the  previous  calcu¬ 
lations  of  the  nonlinear  evolution  of  the  most  unstable  Gbrtler  mode.  The  work  of 
Denier  &  Hall  (1993)  demonstrates  that  the  evolution  of  this  mode  results  in  a  self 
induced  flow  reversal  at  some  finite  distance  downstream  of  the  position  at  which  a 
disturbance  is  first  introduced  into  the  boundary  layer.  This  result  was  independently 
confirmed  by  Otto  &  Bassom  (1993)  in  their  work  on  the  effect  of  crossflow  on  the 
nonlinear  evolution  of  the  most  dangerous  Gbrtler  mode.  However,  in  both  cases,  it 


was  demonstrated  that  the  total  flow  field  becomes  highly  inflectional  in  nature  prior 
to  the  point  at  which  flow  reversal  is  found. 

To  consider  the  secondary  instability  of  such  highly  inflectional  profiles  we  consider 
perturbations  to  the  total  velocity  field  in  the  form  of  inviscid  Rayleigh  waves.  The 
spatial  and  temporal  scales  of  these  modes  are  on  the  same  order  as  the  viscous  wall 
sublayer  in  which  the  nonlinear  vortex  resides,  namely  0(G~}/5  Re~1'2).  We  consider 
perturbations  to  the  nonlinear  vortex  state  of  the  form 

A  ( 'G*U,G*V,G*W,P )  e^«1/5^1/J(x-cO? 

where  tv  is  the  streamwise  wave  number  and  c  is  the  complex  wave  speed;  here 
U  ,V  ,W ,  F  are  functions  of  the  scaled  variables  y  and  z  and  A  is  the  small  per¬ 
turbation  parameter.  The  equation  governing  the  perturbation  quantities  are  most 
easily  written  by  eliminating  C7,  V,  W  to  give  the  three-dimensional  Rayleigh  pressure 
equation 

d2P  d2P  2  - _ 2 

dy2  ^  dz2  01  u  — 

We  impose  the  usual  inviscid  boundary  conditions,  namely  vanishing  normal  velocity 
at  the  solid  boundary  together  with  the  requirement  that  the  perturbation  decays  as  we 
leave  the  viscous  wall  layer.  In  terms  of  the  pressure  perturbation  these  requirements 
become 

dP 

—  =  0  at  y  =  0,  P  — >  0  as  y  —*  oo.  (2.2) 

dy 

We  also  impose  the  condition  that  the  wave  has  the  same  spanwise  period  as  the 
underlying  nonlinear  vortex  state,  hence 

P{z)  =  P(z  -f  2tt),  (2.3) 

where  k  is  the  spanwise  wavenumber  of  the  vortex  velocity  field.  The  form  of  0  is 
given  by  the  total  streamwise  velocity  from  the  vortex  calculation,  and  is 

u(x,y,z)  =  y  +  U(x,y,z),  (2.4) 

with  u(z)  —  u(z  +  27t);  the  term  y  in  relation  (2.4)  is  the  basic  shear  in  the  viscous 
wall  layer  and  U  is  the  vortex  induced  flow.  This  system  now  represents  an  eigenvalue 
problem  which  can  be  solved  for  a  given  a  to  determine  a  complex  phase  speed  c.  and 
thus  determine  the  inviscid  temporal  instability  of  the  flow  at  a  given  streamwise,  x 
location. 


du  dP  du  dP 
dy  dy  ^  dz  dz 


3.  Numerical  techniques 


In  this  section  we  describe  the  techniques  used  to  solve  the  elliptic  equation  (2.1 ). 
We  discretize  the  system  using  a  five-point  regular  spanwise  stencil  and  a  stretched 
three  point  stencil  in  the  normal  coordinate.  As  the  Rayleigh  waves  are  presumed  to 
have  the  same  period  (or  integer  divisors  thereof)  as  the  vortex,  we  alias  the  point 
at  z  —  2ir  to  have  the  same  value  as  that  at  z  =  0.  The  stretching  in  the  normal 
coordinate  ensures  that  the  far  field  boundary  condition  can  Ire  imposed  at  a  suitably 
large  (but  finite)  value  while  still  retaining  resolution  at  the  wall.  The  discretized 
version  of  equation  (2.1)  may  be  written  as 

Ai'Pj+i  +  Bi’Pi  +  C  =  0,  (3.1) 

where  V,  ~  [Pi\,Pii,  ■  ■  •  PiM^  >  where  the  subscript  i  denotes  the  value  at  y,.  In 
the  system  (3.1)  the  matrices  A,  and  C,  are  diagonal,  a  fact  that  is  exploited  by  the 
particular  block  solver  written  to  solve  (3.1).  To  impose  the  asymptotic,  condition  we 
notice  that  P  satisfies 

d*p  2  dP  7fi  n 

- a'P  -  0, 

dy1  y  -  c  dy 

as  y  — +  oo,  where  we  have  made  use  of  the  fact  that  the  vortex  is  confined  to  the 
viscous  wall  layer.  This  equation  has  the  decaying  solution 

P  =  {y  +  A)e~a\ 

where  A  is  a  function  of  a  and  c  but  is  not  relevant  here.  We  choose  to  impose  this 
condition  at  y  —  yoo  by  using  the  Robin  condition, 

1  dP  1 

P  dy  y 

This  outer  limit  is  chosen  so  that  changing  it  does  not  affect  the  calculation.  In  this 
study  we  choose  to  allow  both  odd  arid  even  modes,  and  thus  we  solve  over  the.  whole 
period  rather  than  half  the  range  as  in  Hall  &  Horseman  (1991).  This  allowed  us 
to  have  mixed  modes  which  are  necessary  in  the  three-dimensional  basic  states  cases. 
The  problem  was  normalized  using  the  same  method  as  in  Hall  &  Horseman  (1991) 
by  imposing  the  constraint 

~  =  1  at  y  =  0  Vz.  (3.2) 

dy 
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We  then  iterate  on  the  complex  phase  speed  r  to  drive  the  complex  number, 

_ 1 _ 

/ Real  (P2|j,=0)  +  ilmag  (-P2|j,=o)  dZ 

to  zero,  where  we  use  a  two-dimensional  real  secant  method.  When  the  system  is 
renormalized  the  condition  (3.2)  is  in  fact  Py  —  0,  at  y  —  0. 

In  general  32  points  per  period  were  used  in  the  spanwise  coordinate  whereas  the 
normal  grid  was  made  up  of  90  points,  with  an  infinity  of  yx  —  50.  In  general  the 
same  number  of  points  were  used  for  these  calculations  as  were  employed  in  Otto 
Bassom  (1993)  for  the  vortex  calculations. 


Figure  1.  Growth  rates  ocC{  of  the  primary  and  secondary  modes 


4.  Results  and  Discussion 


We  will  limit  our  discussion  to  the  secondary  instability  of  the  nonlinear  vortex 
state  obtained  by  Denier  Hall  (1993);  further  results  concerning  the  effect  of  crossflow 
on  the  secondary  modes  will  be  presented  in  an  future  paper  (see  Otto  &  Denier  (1993) 
for  full  details).  The  details  of  the  solution  of  the  governing  equations  for  the  nonlinear 
vortex  state  can  be  found  in  Denier  &  Hall  (1993)  and  Otto  &  Bassom  (1993):  the 
reader  is  referred  to  the  aforementioned  papers  for  a  discussion  of  the  numerical  scheme 
used  to  integrate  the  nonlinear  vortex  equations. 

In  figure  1  we  present  the  temporal  growth  rate  of  the  two  major  n.odcs  which  were 
found.  The  modes  marked  ‘primary'  are  conjectured  to  be  projections  of  the  primary 
vortex  onto  the  reduced  inviscid  equations.  The  disturbances  marked  as  ‘secondary' 
are  the  modes  that  we  are  interested  in.  They  have  substantially  greater  growth  rates 
and  thus  are  likely  to  be  more  dangerous.  The  latter  modes  are  characterized  by  the 
amount  of  energy  present  in  the  second  harmonic. 


Pressure  contours  for  a=0.16,  primary 
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Figure  2a.  Iso-pressure  contours  for  the  case  a  —  0.16  (primary). 
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Pressure  contours,  cr  =  0. 162,  secondory 
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Figure  2b.  Iso-pressure  contours  for  the  case  o  =  0.16  (secondary). 

In  figure  2a  we  show  the  contours  of  |P|  for  the  case  o  =  0.16  for  the  primary 
disturbance,  it  should  be  noted  that  this  mode  seems  to  reside  near  the  vortex.  In 
figure  21)  we  show  the  contours  of  pressure  for  the  case  a  =  0.16  for  the  true  secondary 
disturbance.  Notice  that  this  mode  has  maxima  just  ‘outside’  the  vortex,  (as  a  in¬ 
creases  this  phenomena  increases  in  clarity).  In  Figure  3a, b  we  present  the  Fourier 
decomposition  of  the  two  modes,  notice  these  modes  are  comprised  virtually  entirely 
of  cosine  components,  and  hence  we  do  not  show  the  sine  coefficients.  In  Figure  3a 
we  note  that  the  majority  is  confined  to  the  the  fundamental  and  the  two-dimensional 
component  (which  is  not  shown),  however  in  figure  3b  the  mode  labeled  the  ‘secondary’ 
disturbance  can  be  seen  to  have  a  significant  second  harmonic  component.  These  re¬ 
sults  lend  credence  to  the  physical  significance  of  the  ‘primary’  and  ‘secondary'  modes 
pictured  in  Figure  2. 
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Coefficients  of  cosine  tor  primary  a-0.162.  n*  1.2, 3, 4 


Figure  3a.  Fourier  decomposition  of  the  pressure  (primary). 


Coefficients  of  cosine  for  secondo;/  0*0.162,  n*1,2,3.4 


Figure  3b.  Fourier  decomposition  of  the  pressure  (secondary) 


By  imposing  periodicity  in  the  spanwise  coordinate  we  were  aide  to  find  two 
distinct  modes,  here  referred  to  as  the  primary  and  secondary  mode  but  were  unable 
to  demonstrate  the  existence  of  any  ‘odd’  modes  which  lie  rr/2  out- of- -phase  to  the 
underlying  vortex  velocity  field.  However,  by  removing  this  requirement  and  instead 
imposing  boundary  conditions  appropriate  to  such  odd  disturbances  we  are  indeed 
able  to  demonstrate  the  existence  of  such  modes;  the  full  results  of  this  study  will  be 
presented  in  a  forthcoming  paper  (see  Otto  &  Denier  (1993)). 

5.  Conclusions 

We  have  shown  that  flow  situations  involving  a  nonlinear  vortex  are  susceptible  to 
inviscid  disturbances  with  large  growth  rates.  One  of  the  major  points  that  should  be 
noted  from  this  work  is  that  the  most  unstable  .econdary  modes  have  significant  second 
harmonic  components,  which  perhaps  should  be  incorporated  in  the  analytic  approach 
to  this  class  of  problems.  A  similar  phenomena  has  been  observed  in  the  secondary 
instability  of  rotating  disk  flows  by  Balachanclar.  Streett  &  Malik  (1992),  the  reader 
is  referred  to  that  paper  for  further  insight  into  the  nature  of  these  disturbances. 

It  is  not  clear  whether  flows  which  are  inflectional  in  nature  in  the  absence  of 
the  vortex  will  be  stabilized  or  destabilized  by  the  introduction  of  a  vortex  state; 
such  a  question  deserves  further  consideration.  Finally,  it  would  be  an  interesting  and 
significant  problem  to  consider  the  effect  of  increasing  the  amplitude  of  the  inviscid 
wave  to  determine  whether  a  vortex/wave  interaction  could  ensue. 
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